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ABSTRACT 

Large excesses of ^"^Ca in certain presolar graphite and siUcon carbide grains give 
strong evidence for ^^Ti production in supernovae. Furthermore, recent detection of 
the ^^Ti 7-line from the Cas A SNR by CGRO/COMPTEL shows that radioactive ^''Ti 
is produced in supernovae. These make the ^'^Ti abundance an observable diagnostic 
of supernovae. 

Through use of a nuclear reaction network, we have systematically varied reaction 
rates and groups of reaction rates to experimentally identify those that govern ^^Ti 
abundance in core-collapse supernova nucleosynthesis. We survey the nuclear-rate 
dependence by repeated calculations of the identical adiabatic expansion, with peak 
temperature and density chosen to be 5.5xlO^K and 10^ g cm~^, respectively, to 
approximate the conditions in detailed supernova models. We find that, for equal 
total numbers of neutrons and protons (r/=0), ^■^Ti production is most sensitive to 
the following reaction rates: ^^Ti(a,p)^^V, a{2a,-iY'^C, ^^Ti(a,7)^^Cr, ^^\{^,-fY^Ci. 
We tabulate the most sensitive reactions in order of their importance to the ^^Ti 
production near the standard values of currently accepted cross-sections, at both 
reduced reaction rate (0.01 x) and at increased reaction rate (100 x) relative to their 
standard values. Although most reactions retain their importance for > 0, that of 
^^V(p,7)^^Cr drops rapidly for rj >0.0004. Other reactions assume greater significance 
at greater neutron excess: ^"^Cia^-iY^O, ^°Ca(a,7)44Ti, 27Al(a,n)30p, ^^^{{a^nf^S. 
Because many of these rates are unknown experimentally, our results suggest the most 
important targets for future cross section measurements governing the value of this 
observable abundance. 
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1. INTRODUCTION 



Radioactive ^"^Ti, produced in core collapse supernovae, is an isotope of extraordinary 
astrophysical significance. Its primary observable effects appear to be three in number and will 
likely grow in coming years. The first is that the relatively large abundance of ^'^Ca-it is the 
second most abundant calcium isotope and the 44th most abundant species overall in solar system 
material-is overwhelmingly due to its synthesis as '^^Ti parent (Bodansky, Clayton, &: Fowler 1968; 



Woosley, Arnett, & Clayton 1973; Timmes et al. 1996 ). The second is that resulting gamma rays 
from young core collapse supernovae are expected to be visible from several Galactic remnants 
(Clayton 1971; Hartmann et al. 1993| ; Leising & Share 1994; Dupraz et al. 1997). '^^Ti decays to 
"^^Sc emitting 67.9 keV (100%) and 78.4 keV (98%) lines. ^''Sc then decays into "^^Ca, which emits 
a 1.157 MeV (100%) line. The search for Calactic ^^Ti gamma-ray lines has been carried out 
by both the gamma-ray spectroscopy experiments on HEAD 3 ( [Mahoney et al. 1992 ), the SMM 
satellite ( [Leising fc Share 1994| ), and by CGRO surveys with COMPTEL ( Pupraz et al. 1997| ). 
Only from the direction of Cas A SNR has the 1.157 MeV "^^Ti gamma line flux been detected by 
COMPTEL ( llyudin et al. 1994 P^he et al. 1995|) , with a flux initially reported at (7.0±1.7) x 10^^ 
7 cm~^ s~^. Subsequently lyudin et al. (1997) have reported CGRO cycle 1-5 results from a 
2x10® s exposure of Cas A. The result is a flux of (4.8ib0.9) x 10~^ 7 cm~^ s^^ and significance 
of detection of 6a. This measurement is reasonably consistent with the CGRO/OSSE three-line 
flux measurements (The et al. 1996) of (1.8ibl.5) x 10~^ 7 cm~^ s~^ and with the preliminary 



result of RXTE/HEXTE AOl & A02 observations ( [Rothschild et al. 1998|) of (1.3±1.2)xl0-^ 7 
cm~^ s~^, at 4% confidence level with the most probable flux of the combined three instrument 



7 cm s 



This flux translates into (1.3, 7.8) x 10"^ Mq of 



measurements being 3.2 xlO 
^■^Ti for a ^^Ti half-life of (66.6, 39.0) yrs respectively ( [^.IburgeF fc Harbottle 199q ; [Meissner[ 
1996[ ), distance of 3.4 kpc to Cas A ( [Reed et ainMEj) and age of 315 yrs ( [Ashworth 1980]) . 



Because the ^^Ti yield probes the dynamics of core collapse supernova nucleosynthesis, and in 
particular, the location of the mass cut, the pre-supernova composition inside ~2 Mq, and the 
maximum temperature and density reached during the passage of the shock wave in the ejecta, 
this detection has generated great enthusiasm for SNR distance measurements ( [Reed et al. 1995 ), 
for 7-ray line observations (lyudin et al. 1997; Rothschild et al. 1998; The et al. 1996[ ), for new 
nuclear laboratory experiments to accurately determine the '^^Ti half-life (Alburger & Harbottle 



1990; Meissner 1996; Norman et al. 1997[ ; Gorres et al. 1998), and for theoretical nucleosynthesis 



calculations in the context of core-collapse supernovae (Woosley & Hoffman 1991; Timmes et al. 
1996t [Nagataki et al. 1997[). 



The third effect is that ^^Ca-enriched silicon-carbide particles extracted from meteorites have 
been identified (Amari et al. 1992; Niftier et al. 1996; Hoppe et al. 1996| ) as presolar particles that 
condensed within supernova ejecta during their first few years of expansion, while ^*^Ti was still at 
its initial value ( Clayton 1975[ ). These grains are of enormous value in probing the dynamics and 



make up of supernova ejecta ( Clayton, Amari, &: Zinner 1997 ). 

In the near future, we can expect very accurate ^■^Ti line flux measurements, precise distances 
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to several young SNRs, and a accurate ^^Ti half-life (recent measurements of ^'^Ti lifetime by 
Norman et al. 1997 and Gorres et al. 1995 seem to converge at a value of (87.7ibl.7) y). This 
reduction of current uncertainties will allow more meaningful comparisons of the ^^Ti mass in 
Cas A and other young SNRs (in particular, Tycho and Kepler) to supernova models. However, 
a remaining factor limits such comparisons between observations of ^*^Ti mass and supernova 
models, namely, the uncertainty in the nuclear cross sections governing the synthesis of '^^Ti. A 
large number of nuclear reactions play a role in ^^Ti production, and the cross sections for many 
of these are only estimated from nuclear models. While these estimated cross sections should 
be fairly accurate in many cases ( Rauscher, Thielemann, &: Kratz 1997 ), there is no guarantee 
that they are so, and only laboratory measurements will provide us with confidence in the yield 
predictions from nucleosynthesis calculations. Cross section values enter linearly into the stellar 
reaction rate. For this reason, and because of ^^Ti's astrophysical significance, we have surveyed 
the relative importance of the nuclear reaction rates that govern its abundance. Our purpose is to 
identify for the nuclear physics community those reaction rates that have the greatest astrophysical 
importance to the ^^Ti abundance and which, consequently, would be appropriate targets for 
future cross section measurements. Using a nuclear reaction network with 378 isotopes, we have 
systematically varied specific reaction rates and groups of reaction rates in an experimental test 
for astrophysical significance. We identify the most important nuclear rates and explain our 
understanding of the overall rate dependencies we have found. Preliminary results were presented 
in Jin et al. (1997). In the following we will refer to the terms "cross section" and "reaction 
rate" somewhat interchangeably. It should be understood that "reaction rates" are calculated 
from either theoretically predicted and/or experimentally measured cross sections. Related to this 
survey is work done by Woosley & Hoffman (1991) in which they varied the peak temperatures, 
densities, and neutron excesses of explosive silicon burning in alpha-rich freezeout to find the 
range of ^''Co and ^^Ti production with standard reaction rates. 



2. Overall Strategy 

2.1. Laboratory Nuclear Astrophysics 

The present work is part of a long-range effort to determine the identity of the nuclear 
reactions that are of greatest importance to nucleosynthesis in stars. The goal is to improve 
the accuracy of nucleosynthesis calculations by encouraging experimental study of key reactions. 
Although this has always been an important part of the science, as embodied for example in 
the Nobel-Prize winning work of W. A. Fowler, it has not always been emphasized with useful 
clarity. Most nucleosynthesis calculations have been concerned with testing the fundamental 
paradigm: does the evolution and explosion of stars reproduce the known elemental and isotopic 
abundances? Most research efforts address more explicitly the astrophysical models than the 
nuclear data base. That data base is usually taken as a given, and so most calculations have used 
the tabulated values as given in order to calculate the nuclear production within specific stellar 
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models. All have been aware in this that the nuclear data base is uncertain, and conclusions have 
been drawn that respect that uncertainty. But almost no nucleosynthesis calculations display 
the dependence of the abundances on the values of the reaction rates used. The community of 
laboratory nuclear physicists devoted to this problem has continued the work of improving that 
data base. Indeed the level of interest in this has been increasing in the past few years, in part 
because of the construction of several new radioactive-ion-beam facilities ( [Lubkin 1997 ), each of 
which has realized that those facilities offer new opportunities to clarify cross sections near the 
Z=N line in the nuclide chart where several important astrophysical processes occur. 

What is often missing is a clear idea of the relative importance of specific reactions for the 
overall problem. Laboratory scientists quite understandably do not wish to pursue all of the large 
numbers of reactions that naturally occur. It has always been the demonstrated importance of the 
exact value for a cross section that has inspired improved measurements of it. This has usually 
been rather clear in those basic cycles instrumental for nuclear power in stars, where the identity 
of the important reactions is visible. The important reactions during advanced burning processes 
responsible for many of the intermediate mass nuclei have not been so closely scrutinized, however. 
During these processes a large number of reactions occur simultaneously, often presenting parallel 
paths to the same result, so that the significance of the value of each is obscured. The size of a cross 
section is no sure guide to the importance of the value attached to it. In certain circumstances 
during advanced burning a huge reaction rate is often opposed by a huge and nearly equal flow 
due to the inverse reaction, nullifying importance for the exact value of that reaction rate for 
that epoch of that process. So complex are those networks that unambiguous responsibility for 
final abundances has been difficult to assign. This is especially true for nuclear processes taking 
place in or near equilibrium. Our approach to clarifying this is to explicitly change the value of a 
specific reaction rate (or rates) and to redo the calculation to test the sensitivity of a final answer 
to its value. We have developed a logical sequence to do this in a practical manner. 



2.2. Four Requirements for Meaningful Measurements 

Because the goal is to ultimately influence measurements of key cross sections, it is useful to 
list explicitly the four requirements that are, at a minimum, needed for this goal to be achieved: 

1. An appropriate astrophysical model of events significant for nucleosynthesis; 

2. An observable from that process, usually an abundance result that is either known or 
measurable; 

3. Dependency of the value of the observable on the value of a nuclear cross section; 

4. Experimental strategy for measuring that cross section, or at least of using measurable data 
to better calculate it. 
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It is important to note in item 1 that the model must be "appropriate" , not necessarily "correct" . 
We know that models of supernovae, for example, are never "correct" . The unknowns from all of 
physics are many. But the model may be appropriate for testing the sensitivity of the observable 
to the value of the cross section in question. The cross-section dependence may be determinable 
with more accuracy than the absolute value of the observable, because the absolute value depends 
upon the model's realism. Indeed, the model's realism will later be tested by the nucleosynthesis 
results it produces once the important reaction rates are known. We stress this distinction that 
has confused many who ask, "Why is it so important to pin down the key cross sections when the 
astrophysical model introduces more uncertainty than the cross sections do?" 

The observable will usually be some abundance in the natural world. It may be the bulk 
abundance in solar system matter. It may be the abundance outside a supernova that has 
just created new nuclei. It may be an isotopic anomaly within presolar grains having isotopic 
composition distinct from solar isotopes. It could be abundances within cosmic rays. Other 
examples exist. The point is that the observable should on good independent grounds be believed 
to measure the process modeled. Combination of 1 and 2 focuses on the cross section dependence 
of the value of an observable as produced in the process 1. It is then a separate question as to 
whether that cross section has a different importance for another observable, or within a different 
process. Such distinctions are necessary and valid. 

The greater the sensitivity of the observable to the value of a specific reaction rate, the more 
important that rate's value is. Not every abundance observable depends sensitively on the nuclear 
reaction rates; for such cases a cross section does not carry significance (for that observable 
and that process). When experiments are difficult and costly it would be discouraging to learn 
afterwards that the effort was without astrophysical significance. Our approach is to explicitly 
vary the reaction rate to measure the dependence of the final abundance of the observable to the 
value of that rate. By that procedure a quantitative measure of its significance can be presented. 

Item 4, the strategy for the measurement, will be the issue decided by experimental teams 
motivated to nuclear research for purposes of nuclear astrophysics. 



3. The Abundance in the Alpha-Rich Freezeout 



The sensitivity of "^^Ti synthesis to variations in unmeasured cross sections must be determined 
for the process that produces it. Point 1 of §2.2 requires an appropriate model for that process, 
which is now known to be the alpha-rich freezeout of matter that was initially in nuclear 
statistical equilibrium or quasi-statistical equilibrium but that freezes out with excess alpha 
particles ( Woosley et al. 1973 ; Thielemann, Hashimoto, fc Nomoto 1990| ; pi'hielemann, Nomoto, fc 
Hashimoto 1996 ). Although ^■^Ti is also produced during silicon burning ( [Bodansky et al. 1968 ) 
and within explosive helium burning as it may occur within He caps atop Type la supernovae 
(Livne & Arnett 1995; Woosley, Taam, fc Weaver 1986| ), it is believed that the alpha-rich freezeout 
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is responsible for most of natural ^^Ca and it is certain that it is responsible for the ^'^Ti gamma 
lines detectable from young Type II remnants. We therefore choose in this work to determine 
the sensitivity of the '^^Ti yield within that process. It still remains to select an appropriate 
model for that process. Actual supernova models are both uncertain and complicated. In figure 
1^ we show contours of ^^Ti yields computed with our network code (described in more detail 
below) as functions of initial peak temperature and density in nearly adiabatically expanding 
matter experiencing an alpha-rich freezeout. The figure shows that the ^■^Ti yield variation is less 
than a factor of 10 for quite large ranges of peak temperatures and densities near the reference 
parameters, which have been suggested by hydrodynamic models. Therefore, to survey the nuclear 
reaction rate sensitivity it will suffice to evaluate that rate sensitivity within parameterized 
expansions that are typical of the alpha-rich freezeout history within the Type II core. We do this 
by adiabatic expansions of pure ^^Si matter initially at Tg = K = 5.5 and mass density 

p = 10^ g cm~'^. These are the same conditions studied by Woosley & Hoffman (1991) for the 
alpha-rich freezeout (see their Table 1). Such conditions imply a ratio of the number density 
of photons to that of nucleons of 0.57, corresponding to a specific entropy s/k 5, which is 
suitable for shock decomposed silicon in the Type II ejecta. We take the initial density to decline 
exponentially with an e-fold time 0.14 s, a typical hydrodynamical timescale, and we assume that 
the photon-to-nucleon ratio stays fixed at 0.57 throughout the expansion; thus, p oc T^. This one 
zone calculation is then allowed to cool until charged particle reactions freeze out (Tg '^0.25). 
Such a schematic calculation, though inadequate to define the actual '^^Ti yield from core collapse 
supernovae, is adequate as the basis for a survey of cross-section sensitivity for that yield. 

One additional parameter set by the presupernova evolution is the fractional neutron excess 
rj of the matter undergoing the alpha-rich freezeout. Roughly speaking its value, which depends 
upon how far thermonuclear evolution has progressed when the shock wave strikes, is r/ ~0 in the 
He core, r/ ~0.002 (Z^/Zq) in the CO core, and rj ~0.006 in the Si core. We first survey at 77=0 in 
section 4. 



The reaction network used for the calculations is that described in Meyer (1995) and Meyer, 



Krishnan, & Clayton (1996). Tabled shows the nuclides included in the network. Within explosive 
silicon burning that network includes many experimentally unknown reaction cross sections. For 
these reactions, the rates are calculated from the SMOKER code ( p?hielemann, Arnould, &: Truraii 



19871) . 



Figure ^ shows the evolution of mass fractions of key species in the "standard" calculation. 
The material begins as pure ^^Si, but quickly breaks down into light particles which then establish 
a large quasi-equilibrium (QSE) cluster ranging from silicon to beyond nickel. For these particular 
conditions, the system does not attain complete nuclear statistical equilibrium (NSE). This is due 
to the slowness of the triple-a reaction in maintaining the abundance (per nucleon) of heavy 
nuclei {A > 12) at the level demanded by NSE. This slowness imposes an extra constraint on the 
equilibrium, which is thus a QSE, not an NSE ( [Meyer, Krishnan, fc Clayton 19*97 ). Within the 



QSE cluster the rates of reactions involving p, n, a, and 7 reactions proceed at the same rate 
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as their inverses, just as in NSE; but in QSE the total number of nuclei in the cluster is not the 
number they would have in NSE. 

As the temperature falls, the ^^Ti abundance drops. This is a consequence of the fact that the 
QSE cluster shifts upward in mass owing to falling temperature in the face of a deficit of heavy 
nuclei relative to NSE demands ( [Meyer, Krishnan, fc Clayton 1996| ). ^^Ti breaks out of the large 
QSE cluster around Tg = 4.3, and its abundance starts to grow as alpha particles in this alpha-rich 
freezeout reassemble into ^^C. Subsequent alpha captures then carry the newly assembled ^^C 
up to higher-mass nuclei, including ^'^Ti. The final frozen ^^Ti mass fraction ends up roughly 
three orders of magnitude greater than its lowest value at the moment it broke out of the large 
QSE cluster. MPEG movies illustrating the nuclear dynamics of the "standard" calculation are 



available for viewing on the world-wide web at http: / /photon. phys.clemson.edu/movies. html 



4. The Reaction-Rate Survey at r]=0 

In order to determine the sensitivity of the ^^Ti yield to specific reaction rates, we carried out 
a systematic survey. The central plan of the survey is to increase or decrease a reaction rate or set 
of reaction rates by a large factor and then to rerun the alpha-rich freezeout calculation described 
above. The associated errors in the nuclear cross sections may be either a uniform increase or 
decrease of their values as a function of energy within the statistical model, or they may be specific 
energy-dependent contributions to the reaction rate. We then compared the resulting ^^Ti yield 
to that in the "standard" calculation in which no reaction rate was changed. A significant yield 
difference indicates a sensitivity to the varied rates. It is important to note that the network code 
automatically computes reverse reaction rates from the (largely known) nuclear masses, forward 
reaction rates, and detailed balance; thus, for example, if the rate for the reaction ^^Cr(a, 7)^^Fe 
is increased by a factor of 100, the code automatically increased the rate for the reverse reaction 
^■^Fe(7, a)^^Cr by a factor of 100. This is in compliance with detailed balance and is necessary to 
allow the network to relax to the appropriate equilibrium when relevant. 

We began with a survey over element number Z. This step multiplies all charge-increasing 
and mass-increasing reaction rates on all isotopes of a single element Z by a factor of 0.01 and ran 
the alpha-rich freezeout. This set of reaction consists of (a,7), (o:,p), (a,n), (p,7), (p,n), (n,7). 
Elements from carbon (Z = 6) to bromine {Z = 35) are tested in this way. We then repeated the 
procedure, but this time using a factor of 100. The results are shown in figure |3[ From this figure 
it is evident that the ^^Ti yield is sensitive to reaction rates on a number of elements, particularly 
on Ti itself and on V. What is perhaps more striking, however, especially in the lower panel, is 
the large number of elements to which the ^^Ti yield is not sensitive. For example, the ^^Ti yield 
(the astrophysically relevant observable) is not sensitive to large (factor of 100) uncertainties in 
reactions on any isotope of carbon through phosphorus! Such a negative result insofar as rate 
sensitivity is concerned is as valuable as a positive result because it increases the robustness of 
calculations of nucleosynthesis in alpha-rich freezeouts. The uncertainties in these cross sections 
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do not jeopardize the calculated ^'^Ti yields. Understanding this prevents incorrect astrophysical 
motivations for nuclear experiments on reactions irrelevant to the observable. 

There is a simple reason for the insensitivity of the ^^Ti yield to reaction-rate variations on 
so many elements. As discussed in the previous section, the ^^Ti abundance builds up as new 
heavy nuclei assemble late in the expansion. A nearly steady nuclear flow, whose magnitude is 
governed by the triple-a rate, carries these new nuclei up to higher mass. By steady flow, we 
mean that the nuclear flow into a nuclear species by capture from lower-mass nuclei equals the 
capture flow out of that species; thus, the abundance of the nuclear species in question does not 
change except over long time. Because the flow is given by a rate times an abundance, an increase 
in the rate is compensated by a decrease in the abundance so that the net nuclear flow upward 
remains unchanged. By analogy, if a smoothly flowing river widens at some point, its depth or its 
speed decreases in order to keep the same number of gallons per second moving downstream. As a 
particular nuclear example, a factor of 100 decrease in the ^^Ar abundance compensates the factor 
of 100 increase in the ^^Ar(a, 7)^^Ca reaction rate in a steady flow. A large uncertainty in the 
^^Ar(a, 7)'^''Ca rate would thus certainly be important for the yield of ^^Ar, but it is essentially 
irrelevant for our chosen observable, the ^^Ti yield. These considerations illustrate the importance 
of clearly defining the observable when discussing the importance of particular nuclear reactions. 

Table § ranks by element the ^^Ti yield variation resulting from the survey over each element. 
The percent change is the increase in ^^Ti yield divided by the reference value, in parts per 
hundred. A large variation for a given element means that the final ^"^Ti yield is sensitive to at 
least one reaction rate on some isotope of that element. The next task is to identify the key 
isotope. For each element showing a large ^^Ti-yield variation, we performed an isotope search by 
multiplying all reactions on a given isotope of that element Z by a factor of 0.01 and rerunning the 
alpha-rich freezeout. This was repeated for each isotope of the element Z in the network. After 
completing the isotope survey for element we surveyed the next important element. This entire 
procedure was then repeated but with a factor 100 multiplying the reaction rates. 

Figure ^ shows the results for the isotope surveys for Ar, Ca, Ti, and V for both the 0.01 and 
100 factors. It demonstrates that the ^"^Ti yield significantly varies for changes on only one isotope 
of each of the elements. For the elements shown, these are ^^Ar, ^''Ca, ^^Ti, and ^^V. Table ^ ranks 
the particular isotopes for the 0.01 and 100 reaction-rate variation surveys. Comparing the first 
entries of tables |2| and |3|, for example, shows that a 100-fold decrease of the reaction rates for all 
Ti isotopes produces an almost identical increase (+373%) in final ^^Ti abundance, in comparison 
to the increase (+372%) for the ^^Ti reaction rate only. 

With the key isotopes now identified, we next determined which particular nuclear reaction 
dominated the sensitivity of the ^^Ti yield. This step of the survey, similar to that done above, 
varied a specific reaction on the key isotope. The first such survey was of the nuclear reactions on 
^^V. We multiplied the ^^V(n, 7)^^V rate by a factor of 0.01 and reran the alpha-rich freezeout. 
We then repeated this for the (p, 7), {p,n), (0,7), {a,n), and {a,p) reactions. We did this for 
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each isotope hsted in table ^. We then repeated the procedure with a multiplication factor 100. 

The results for the surveys on ^'^Ca, ^^Ti, ^^Ti, and ^^Fe are shown in figure ^. Its discrete 
abscissa is the reaction channel for that isotope. The results clearly show that the ^^Ti yield 
is sensitive to variations in the rates of the particular reactions ^'^Ca(a, 7)'^^Ti, '^^Ti(a,p)^'^V, 
^^Ti(a, 7)^^Cr, ^^V(p, 7)^^Cr, and ^^Fe(a, n)^'^Ni. Table |^ ranks the particular nuclear reactions 
according to their influence at r]=0 on the final ^^Ti yield. For example, reduction of ^^Ti(a, 7)^^Cr 
has almost no effect (unless of course the {a,p) branch is also reduced), whereas increase of its 
rate by 100 x is almost as effective as an increase of the {a,p) channel. Such conclusions depend 
on the relative magnitude of proton and gamma widths of the compound nuclear states, and on 
the factor by which they are altered. Owing to its special role in the late assembly of heavy nuclei, 
we also surveyed the triple-a reaction and include it in table ^. 

One might have predicted that a number of these reactions would be important, although the 
magnitude of their effects would not have been as easy to guess. The ''5V(p,7)^6Cr reaction, the 
one showing the largest effect, however, was not one for which we had anticipated a special role. 
This underscores the importance of making a systematic survey when considering the reaction-rate 
sensitivity of a particular observable. In fact, not only are systematic surveys important for clearly 
delineating reaction-rate sensitivities, but they can also clarify nucleosynthesis. It was only after 
developing new insights into the nuclear dynamics of an alpha-rich freezeout that we were able to 
understand the reason for the importance of the ^^V(p,7)^^Cr reaction. By the same token we 
can understand why its importance vanishes for neutron richness rj >0.0004 (see Section 6). We 
return to this issue in the next section. 

The remaining task for the ^^Ti survey was to determine the form of the sensitivity of each 
particular important reaction. Figure ^ shows that increasing the ■*^V(p, 7)^^Cr reaction by a 
factor of 100 decreases the resulting *^^Ti yield by a factor of nearly 50. This suggests a nearly 
linear dependence of the ■^^Ti yield on the rate of this reaction. For a decrease by a factor of 100 
on this rate, however, the "^^Ti yield increases by less than a factor of two, suggesting a non-linear 
dependence. The yield dependence on the rate value determines its significance. To quantify 
the form of the dependence, the alpha-rich freezeout was repeated many times, each time with 
a specific important nuclear reaction rate multiplied by a factor between 0.01 and 100. This 
maps out the sensitivity of the ^■^Ti yield to any change by a factor between 0.01 and 100 on the 
important nuclear reactions. Figures |6|, 0, and ^ show results of such surveys. In each of these 
figures, a small solid square indicates the result using the "standard" value for that particular rate. 

These figures illuminate more subtle aspects of the reaction rate sensitivities. Figure |6| shows 
that a factor of two increase or decrease in the ^^V(p, 7)^^Cr rate results in a ~ 20% change in 
the yield of ^"^Ti. On the other hand, figure ^ shows that a factor of two change in the ^^Ti(a, 7) 
reaction rate barely changes the ^^Ti yield, even though a factor of 100 increase in the rate 
dramatically drops the yield. Because the reaction rates we use are either experimentally known 
or among the best theoretical estimates presently available, we expect the true rates in general 
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not to be more than, say, an order of magnitude different. Thus, the quantity expressing the 
^^Ti-sensitivity to a particular reaction near its reference value is the slope of the curve for that 
reaction, in analogy to figure ^. In table ^ we present these quantities, in a ranking by magnitude, 
for the sixteen most important reactions at r/=0 insofar as final ^'^Ti abundance is concerned. 

As an example of how to use table ^, consider the case of ^^V(p, 7)^^Cr. If it were measured, 
one would then compare that rate (especially near Tg = 2 for reasons to be discussed in the next 
section) to the value used in our surveys, the fits for which we present in table |8[ If the measured 
rate were a factor of 1.5 greater than the tabulated value for all temperatures, the resulting ^^Ti 
yield would change by -0.361 x (1.5 - 1) = -0.18 = -18%, that is, it would be 82% of the 
standard value. 

It was only through a large and systematic survey that we were able to determine the 
relative importance of the reactions listed in table ^. This was a computationally intensive effort 
comprising more than 2000 sets of alpha-rich freezeout calculations, each taking ~20 minutes on 
a 200 MHz DEC Alpha workstation. Nevertheless, it is only in this way that one can be fully 
confident of understanding the reaction-rate sensitivity of the ^"^Ti yield. 

A final part of the survey applies to neutron-excess tj greater than zero. The survey is repeated 
in summary form in Section 6 for rj = 0.002 and rj = 0.006. For these neutron enrichments, ^^Ti 
synthesis was unaffected by the value of the '^^V(p,7)^^Cr reaction rate, a result also explained in 
Section 6. But first we describe the reasons for the importance of the Table ^ reactions at r/ = 0. 



5. Discussion of the Most Important Reactions 

Table ^ lists the reactions to which the ^^Ti yield is most sensitive in an alpha-rich freezeout. 
In this section, we briefiy explore the reasons for this sensitivity for three important reactions. 
In addition to the immediate interest of understanding the importance of these reactions, this 
exercise provides new insights into the nuclear dynamics of the alpha-rich freezeout. 

Any nucleosynthetic system seeks to maximize its entropy and thereby to achieve an 
equilibrium. Its ability to do this, however, is governed by the rate of certain key reactions. When 
a particular reaction becomes too slow, the equilibrium breaks. The system still maximizes the 
entropy, but the slowness of the particular reaction adds a constraint which restricts the number 
of accessible states (Meyer, Krishnan, & Clayton 1997). As the system evolves toward lower 



temperatures, more reactions become slow, more equilibria break, and more constraints are added 
to the entropy maximum. In this way, the system descends the "hierarchy of statistical equilibria" 
(contribution of Meyer in Wallerstein et al. 1997 ). Several of the reactions in table |^ are important 



because they control the breaking of equilibria and, thus, the descent of the hierarchy. Other 
reactions are important because they control when the steady (but non-equilibrium) flow into ^^Ti 
breaks down or because they control the economy of light particles that influence the abundance 
of ^^Ti late in the expansion. 
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5.1. The triple-a reaction 

We begin with the triple-a reaction. Figure IC shows the time history of the *^^Ti mass 



fraction during two alpha-rich freezeouts-one with the standard ( [Caughlan fc Fowler 1988D rate 
(dashed curve) and one with the standard rate increased by a factor of 100 (solid curve). For 
comparison, the dotted curve gives the ^^Ti mass fraction in NSE. This NSE is calculated at 
each timestep in the expansion using the same temperature, density, and neutron richness as in 
the network. The dashed and solid curves differ throughout the expansions. In the case with the 
increased triple-a rate, the nuclear system achieves and maintains NSE for the early part of the 
expansion. With the standard triple-a rate, the system never reaches NSE but rather only a QSE 
in which the number of heavy nuclei Yh differs from that in NSE. Below Tg ~ 5.3 the number of 
heavy nuclei in the standard run QSE is less than that in NSE, so the heavy nuclei are in the 
presence of an overabundance (relative to NSE) of light particles. This favors nuclei heavier than 
"^^Ti (e.g. ^^Ni), so the ^^Ti mass fraction is far below that of NSE until later in the expansion 
when the large QSE breaks and the ^^Ti can build back up without being driven upward to ^^Ni. 

The ^*^Ti mass fraction in the standard calculation is less than that in the case with the 
increased triple-a reaction for Tg < 5.3. The latter expansion deviates from NSE beginning at 
Tg ~ 5.2. At this point the nuclear system is also in a QSE, but this QSE has a greater number 
of heavy nuclei and a smaller number of light particles than the standard case QSE. This yields 
a greater mass fraction of ^^Ti within the QSE cluster. At the point ^^Ti breaks out of the large 
QSE cluster containing ^^Ni (the dip in both curves near Tg = 4), the ^^Ti mass fraction in the 
expansion with the increased triple-a rate exceeds that in the standard case by a factor of ~ 
eight. More importantly, the expansion with the increased triple-a rate has a factor of ~ 4.6 lower 
mass fraction of ^He. Remarkably, this factor of 4.6 (when cubed) compensates for the increased 
triple-a cross section; thus, the two expansions create new nuclei at the nearly the same rate at 
late times. The crucial aspect for ^^Ti is that the lower "^He abundance in the expansion with the 
increased triple-a cross section requires larger abundances of nuclei between ^^C and ^^Ni to carry 
the same flow. This causes the increase in the final yield of ^'^Ti. 



5.2. 44xi(a,p)4^V 



Figure 11 shows the ^^Ti mass fraction in the standard calculation (dashed curve) and in an 
expansion with the ^^Ti(a,p)'^^V rate increased by a factor of 100. Also shown as the dotted curve 
is the QSE mass fraction of ^^Ti in the standard expansion calculated at each timestep from the 
same temperature, density, neutron richness, and number of heavy nuclei as in the network. The 
evolution of the ^^Ti mass fraction follows QSE in both cases down to Tg ~ 4.5. The increased 
^^Ti(a,p)^'^V rate, however, allows the ^^Ti to remain in QSE with ^^Ni longer; thus, when the 
QSE finally breaks, the ^^Ti mass fraction is ~ 30 times lower than the corresponding value in the 
standard case. This result illustrates the important role ^^Ti(a,p)^'^V plays in linking the Si-Ca 
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QSE cluster to the Ni-centered QSE cluster ( Woosley et al. 1973| ; [Hix fc Thielemann 1996 ). Once 
the two clusters fall out of equilibrium with each other, the ^^Ti(a,p)'^''V reaction governs the rate 
at which abundance moves from the Si-Ca-Ti cluster to the Ni cluster. A faster rate gives a lower 
abundance in the Si-Ca-Ti cluster and a lower ^'^Ti yield. 



5.3. 45v(p,^)46cr 

Perhaps the biggest surprise in our survey was the importance of the '^^V(p, 7)^^Cr reaction. 



Figure |T2| shows the "^"^Ti mass fraction for the standard calculation (dashed curve) and for an 
expansion with the ^^V(p, 7)^^Cr increased by a factor of 100. The ^^Ti evolution is precisely the 
same in the two expansions down to Tg ~ 2. The subsequent deviation between the two cases, 
however, is dramatic. While the ^^Ti mass fraction drops by ~ 30% as Tg falls below 2 in the 
standard calculation, it plummets by a factor of ~ 100 for the increased ^^V(p, 7)^^Cr rate. The 
reason is, as in the previous cases, that this reaction controls the breaking of an equilibrium. 

The equilibrium of importance in this case is a (p, 7) — (7,p) equilibrium among the N = 22 
isotones. In such equilibrium, the rates (as opposed to the cross sections) of (p,7) reactions proceed 
at the same rates as the reverse (7,p) reactions; i.e. "^Z -|- p ^ "^+-^(Z+1) + 7, where the reversed 
arrows indicate equal numbers of reactions per unit time. We focus on the N = 22 isotones, A = 
Z + 22. Figure |l^ shows how this equilibrium breaks for the standard set of cross sections. It 
gives for four temperatures the ratio of the network abundance of these isotones to the values they 
would have in (p, 7) — (7,^) equilibrium. For Tg = 2.66, all of the isotones for Ti through Fe are in 
equilibrium with each other under exchange of protons, but the lower-charge isotones have already 
broken out of this equilibrium. As the temperature falls to Tg = 2.15, '^^Fe strongly falls out of 
equilibrium with the other four species, but also evident is the fact that the two small equilibrium 
clusters (^^Ti & ^^V) and ('^^Cr & ^''Mn) are starting to fall out of equilibrium with each other; 
that is, ^^V and ^^Cr are no longer in (p, 7) — (7,^) equilibrium with each other. This break grows 
with decreasing temperature. By Tg = 1.81 the two small (p, 7) — (7,p) clusters, (*^^Ti & ^^V) and 
(^^Cr & ^''Mn), are strongly out of equilibrium with each other. When this occurs, ^^V(p, 7)^^Cr 
produces net destruction of ^^V, and hence of '^^Ti, with which it is in equilibrium. 

The Tg = 1.81 panel of this figure shows that the '^^Ti-^^V cluster is overabundant and the 
46Qj-_47]y[j^ cluster is underabundant with respect to the overall (p, 7) — (7,p) equilibrium. This 
indicates that too few nuclei have moved from the 4*^Ti-^^V cluster up to the ^^Cr-^'^Mn cluster 
because of the slowness of the reaction that links them, namely, '^^V(p, 7)^^Cr. By increasing the 
rate of this reaction, the two small clusters remain in equilibrium with each other longer, the 
N = 22 isotones shift to higher mass, and more '^^Ti is destroyed. 

The dotted curve in figure |l^ amplifies this point by showing what would happen in the 
standard case to the 4*^Ti mass fraction if the (p, 7) — (7,p) equilibrium among N = 22 isotones 
persisted to low temperature. The faster ^^V(p, 7)^^Cr reaction allows the ^^Ti to remain in the 
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(PiT) ~ (ItP) equilibrium to lower temperature. Because the (p, 7) — (7,p) equilibrium abundances 
shift to higher mass with lower temperature, the ^'^Ti mass fraction falls. Even with the fast 
45v(p, 7)'^^Cr reaction, however, the equilibrium eventually breaks and the ^^Ti mass fraction 
freezes out. In the surveys at greater neutron excess to follow, we demonstrate that the controlling 
significance of ^^V(p,7)^^Cr vanishes for 77 > 0.0004. This will be understood as a consequence of 
the smaller free proton density with increased rj. 



6. Surveys Having Excess Neutrons 

We repeated the surveys in exactly the same procedure for a nuclear gas containing more 
neutrons than protons in total. This initial neutron richness, when the shock wave strikes the 
overlying matter, is not altered during the subsequent rapid nuclear burning, which occurs too 
rapidly for significant electron capture to occur. Thus the gas can be described by a constant 
parameter t] = J^ii^i ~ Zi)Yi, where Ni and Zi are the mass fractions of nucleons (both bound in 
species i and free) and Yi = Xi/Ai is the abundance of species i with mass fractions Xi and atomic 
mass number Ai. The survey shows that although most nuclear reactions retain their importance 
for 4'^Ti production during alpha-rich freezeouts, the ^^V(p,7)^^Cr reaction is very sensitive to the 
value of r]. Demonstrating and understanding these facts is the main goal of this section. 

The range of anticipated values for r/ is set by the presupernova evolution and structure and 
by the strength of the shock wave. The latter determines how far radially the shock can propagate 
with sufficient strength to establish an alpha-rich freezeout, without which the yield of ^"^Ti is 
too small to contribute substantially to the natural ^^Ca abundance. Within the entire He core 
of massive stars the matter is dominated by He, ^"^N, ^^O and ^^C, all of which have r/ = 0. 
But within the CO core r/ has been increased by the conversion of ^^N to ^^O and ^^Ne, both of 
which have two excess neutrons. At solar initial metallicity {Zq) these comprise about 2% by 
mass, so that rj = 0.002 throughout the CO core. Supernovae of lower initial metallicity Zj, such 
as supernova 1987A in the Magellanic Cloud, will have 77 near O.OO2(Zj/Z0) in their CO cores. 
Hydrostatic carbon burning does not much increase r/ within the NeMg core, but hydrostatic O 
burning does do so. The more central core in which O has burned to create a Si core has acquired 



more excess neutrons by electron captures on the products of O burning (Woosley, Arnett, & 



Clayton 1972] ), yielding r/ in the range 0.003-0.006. In the interests of limiting the computational 



demands, we have performed the full survey at the values 77 = 0.002 and 0.006. Finer precision 
mapping in was also done for the ^^V(p,7)^^Cr reaction alone, after determining its sensitivity 
to the neutron excess. 

Tables ^ and list the effect of reaction rate changes for the ten most important 
regulators of ^^Ti production at these two values of ry. Many of the same reactions remain 
important, such as ^^Ti(a,p)^'^V, '^''Ca(a,7)^'^Ti, and '^^Ar(o;,p)^^K; but the importance of 
^^V(p,7)'^^Cr disappears and some new important reactions appear (e.g. ^^C(a,7)^^0). More 
complete lists can be read and downloaded from the Clemson nuclear astrophysics web site 
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( http: / / photon.phys.clemson.edu / tables.html ) . 



In the interest of understanding we comment upon two of these dramatic changes. The 
importance of the ^^V(p,7)^^Cr reaction decUnes precipitously for t] > 0.0001, as shown in figure 



14. The cause is the decline of the free-proton abundance Yp at every temperature with increasing 
7]. As T] increases, protons tend to be more tightly bound in nuclei. Although the ^^Ti-^^V pair 
remain in (p, 7) equilibrium down to Tg^ 2 for r] as large as ~ 0.01, the ^^V abundance is declining 
with increasing r] due to the diminishing free proton abundance. As a result the final proton 



captures are unable to reduce the abundance of the ^^Ti-^^V pair during the freezeout. Figure 14 
shows the free proton abundance at Tg = 2, the temperature at which the (p, 7) equilibrium begins 
its freezeout, for several calculations at different r/'s. In the equilibrium, the abundance ratio 
'^^V/^^Ti oc Yp. The ^^V(p,7)^^Cr reaction rate is also proportional to Yp] thus, the importance of 
this reaction rate in modifying the final ^^Ti yield goes as Yp. The fact that Yp at Tg =2 drops 
off so rapidly for rj > 0.0001 explains why the final ^^Ti yield sensitivity also declines so quickly. 
The consequence is that the ^^V(p,7)'^^Cr reaction loses considerable importance for observations 
of ^^Ti in supernovae, at least unless they have very small initial metallicity, or unless the shock 
wave is able to reach the He shell, or unless some other process keeps the proton abundance higher 
than would be expected simply from the alpha-rich freezeout. It is worth noting that some Type 



la events create abundant ^^Ti ( Clayton et al. 1997 ); but the burning is not literally an alpha-rich 



freezeout although it does have hot excess alpha particles. A variety of r/'s might be expected in 
this scenario depending on the initial composition. For low r] we will again expect certain proton 
capture reactions to affect the ^"^Ti yield. We will study the important cross sections for many 
observables in this setting in a subsequent work in this series. 

The ^^C(a,7)^^0 reaction assumes greater importance for 77 > 0. The greater neutron richness 
of the matter causes alpha particles to bind more tightly into nuclei, thereby decreasing the alpha 
particle abundance throughout most of the expansion. Because of the smaller alpha particle 
abundance, flow of newly assembled nuclei to higher mass freezes out at higher temperature. 
When the cross section for alpha particle capture on ^^C is decreased, this freezeout occurs even 
earlier, leading to considerably reduced ^^Ti production. 



7. Conclusion 

We have analyzed the influence of individual reaction rates on ^^Ti production in alpha-rich 
freezeouts within Type II supernova events. Using direct surveys we have established that the ^^Ti 
production at r/ = is most sensitive to the following reaction rates: ^^V(p,7)^^Cr, a(2a,7)-'^^C, 
44Ti(a,p)47V, ^^Ti{a,-f)^^CT, ^^Ca{a,-f)^^Ti, and 57Co(p,n)57Ni. For 7/ > the two proton-induced 
reactions decline rapidly in importance and some other reactions become crucial: ^^C(a,7)^^0, 
'^°Ca(a,7)^''Ti, ^'^Al(a,n)^°P, and ^°Si(a,n)^^S. Several good reasons justify this effort at a depth 
of detail and computational time that exceeds that normally devoted to nuclear studies in 
astrophysics. Three main reasons are: (1) to provide guidelines for laboratory measurements; 
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(2) to establish better diagnostics of supernova events; (3) to gain insight into comphcated and 
nonhnear nuclear reaction networks of importance for nucleosynthesis. 

The first reason parallels the large increase in recent activity in laboratory nuclear astrophysics. 
In particular, many new radioactive- ion-beam facilities are rapidly establishing the capability 
to measure reaction cross sections near and on the proton-rich side of the Z = N line. This 
capability opens to experimental study nucleosynthesis processes that have heretofore had to rely 
on computed cross sections. Each such experiment will be costly and time consuming, however, 
so that solid grounds for the importance of any specific reaction are welcome to experimental 
planners. 

The second reason parallels the large increase in recent observations of supernova 
nucleosynthesis. Although the mass of ^^Ti within the Cas A supernova 310-year-old remnant 
cannot yet be regarded as measured with good precision, the recent detection of it insures that it 
eventually will be well measured. When that precision has been achieved, the mass of '^^'pj ^jj^ 
reflect the alpha-rich-freezeout mass that was ejected. Since this mass depends in exciting ways on 
the physics that underlies the Type II event (the mass cut and the shock mechanism), all manner 
of exciting issues may hinge upon its production. The same will be true for the several recent 
Galactic supernovae that will be detected with the improved gamma-ray detector sensitivity that 
now exists in laboratories. It is probable that significant variations of the ^^Ti yield exist between 
different events. Once that observational goal is achieved, the most accurate cross sections will 
be wanted so that the scientific uncertainties reflect the uncertainties over the realistic supernova 
model rather than over the nuclear cross sections. Similar excitement attends the measurement of 
44rp-^48'pj production ratios in samples of supernova matter as recorded by supernova condensates 
found in meteorites. One does not want nuclear cross section uncertainties to stand in the way of 
quantitative analysis. 

The third reason parallels the large recent increase in the sophistication with which 
nucleosynthesis is viewed. The present study illuminates this very well. Within the complex 
networks of explosive nucleosynthesis, one can ideally think of the importance of any specific 
nuclear reaction in various ways. It may have almost no significance if that reaction maintains a 
partial QSE throughout the burning. When QSE breaks down or fragments, a reaction may play a 
rather complicated role in governing the changing relative total numbers within the different QSE 
clusters; or it may play a role in the density of free light particles; or it may simply be the specific 
channel by which an abundance is altered in the freezeout. The important reactions identified by 
us for '^^Ti production in the alpha-rich freezeout reveal each of these aspects. We have discussed 
these to some extent in the discussions of the important reactions. 

Finally we woTild observe that our study concentrated on one very important nucleus within 
one nucleosynthesis process. There are many nucleosynthesis processes producing many different 
nuclei, so additional studies of the present type are recommended. This requires judgment as 
well as computing power. Clearly one does not wish to have huge numbers of irrelevant studies 
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performed. For this reason, such studies should follow the "four-requirement" structure outlined 
in §2.2. 

We would like to thank the referee, Robert Hoffman for a thorough review and suggestions to 
improve this paper. This work has been supported by NASA grants NAG5-4329 and NAGW-3277. 
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Table 1. Nuclei in the Network 
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Table 2. Order of importance of elements affecting the Ti production in alpha-rich freezeout 
mechanism if the nuclear cross sections for all mass or charge increasing reaction channels on all 
isotopes of an element are changed by a factor of 1/100 or 100. 



Rank 


Reaction rate 
multiplied by 1/100 
Element % ^^Ti change 


Reaction rate 
multiplied by 100 
Element % ^Ti change 


1 


Ti 


+372.6 


V 


-98.0 


2 


Ca 


-72.6 


Ti 
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Table 3. The Order of Importance of Isotope to the ^^Ti final mass fraction when all mass or 
charge increasing nuclear reactions of the isotope are multiplied by 1/100 or 100. 



Rank 


Reaction rate 


Reaction rate 




multiplied by 1/100 


multiplied by 100 




Isotope 


% ^^Ti change 


Isotope 


% ^^Ti change 


1 




+371.7 


45v 


-98.0 


2 


«Ca 


-72.5 


44Ti 


-91.2 


3 


45y 


+56.6 


s^Co 


+27.1 
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36 Ar 


-49.7 


40Ca 


+22.0 
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-49.5 
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54pg 


+9.4 


7 


32s 


-20.4 


36 Ar 


+5.6 


8 


13N 


-14.8 


46 Cr 


-5.3 


9 


58CU 


-13.5 


56Ni 


-5.2 


10 


56Ni 


-4.5 


61CU 


+4.0 


11 


28Si 


-4.2 


58Ni 


+3.4 


12 


^6 Co 


-4.1 


23Na 


-2.6 


13 


20Ne 


-4.0 


12c 


-1.8 


14 


12c 


+3.1 


13N 


-1.5 


15 


48Cr 


+2.8 
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Table 4. Order of importance of reactions producing ^'*T[ at r/=0 from low multiplication factor 

(O.Olx) and high multiplication factor (lOOx). 



Rank 


Reaction rate 




Reaction rate 




multiplied by 1/100 




multiplied by 100 




Reaction % ^^Ti change 


Reaction % ^Ti change 


1 


44Ti(a,p)47V 


+173 


^^V(p,7)^6Cr -98 


2 


a(2a,7)i2c 


-100 


a(2a,7)i2c +67 


3 


40Ca(a,7)44Ti 


-72 


44Ti(a,p)47v .89 


4 


45v(p,7)46Cr 


+57 


^Ti(a,7)48cr -61 


5 


57Ni(p,7)58Cu 


-47 


57Co(p,n)57Ni +25 


6 


^^Co(p,n)5^Ni 


-33 


40Ca(a,7)44Ti +22 


7 


13n(p,7)140 


-16 


^^Ni(n,7)5SNi +10 


8 


58Cu(p,7)59Zn 


-14 


5%e(a,n)57Ni +9.4 


9 


36Ar(a,p)39K 


-11 


36Ar(a,p)39K +5.5 


10 


i2C(a,7)i60 


+3.5 


36Ar(a,7)40Ca +5.3 
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Table 5. Order of Importance of Reactions Producing ^^Ti at r/ = According to the Slope of 

X(^'^Ti) Near the Standard Cross Section 



Reaction 


Slope 


^^Ti(a,p)4^V 


-0.394 


a(2Q,7)i2C 


+0.386 


«V(p,7)*^Cr 


-0.361 


^°Ca(a,7)^^Ti 


-F0.137 


5''Co(p,n)«^Ni 


+0.102 


3«Ar(a,p)39K 


+0.037 


**Ti(a,7)*«Cr 


-0.024 




-0.017 


^•^Ni(p,7)5^Cu 


+0.013 


5*Cu(p,7)5«Zn 


+0.011 


3«Ar(Q,7)«Ca 


+0.008 


^^Ti(p,7)*^V 


-0.005 


5^Co(p,7)58Ni 


+0.002 


5''Ni(n,7)5«Cu 


+0.002 


5*Fe(a,n)5^Ni 


+0.002 


*°Ca(a,p)«Sc 


-0.002 
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Table 6. Order of importance of reactions producing Ti at ry=0.002 from low multiplication 
factor (0.01 x) and high multiplication factor (lOOx). 



Rank Reaction rate Reaction rate 

multiplied by 1/100 multiplied by 100 

Reaction % ^'^Ti change Reaction % ^^Ti change 



1 


4^Ti(a,p)47v 


+208 


44Ti(a,p)47V 


-93 


2 


i2C(a,7)i60 


-72 


44Ti(a,7)48Cr 


-66 


3 


40Ca(a,7)44Ti 


-66 


27Al(a,n)30p 


-60 


4 


20Ne(a,7)24Mg 


-16 


30Si(a,n)33s 


-33 


5 


30Si(p,7)3ip 


-9.2 


i2C(a,7)i60 


+18 


6 


36Ar(a,p)39K 


-7.9 


40Ca(a,7)44Ti 


+15 


7 


59Ni(p,n)59Cu 


-4.7 


23Na(a,p)26Mg 


-4.7 


8 


5'^Ni(p,7)60Cu 


-4.7 


39K(a,p)42ca 


+4.7 


9 


44Ti(a,7)48cr 


+2.8 


27Al(p,7)28Si 


+4.3 


10 


27Al(a,n)30p 


+2.7 


24Mg(a,7)28Si 


+4.2 



Table 7. Order of importance of reactions producing *^Ti at 77=0.006 from low multiplication 
factor (0.01 x) and high multiplication factor (lOOx). 



Rank Reaction rate Reaction rate 

multiplied by 1/100 multiplied by 100 

Reaction % ^'^Ti change Reaction % ^'^Ti change 



1 




+211 




-93 


2 


i2C(a,7)i60 


-79 


4^Ti(a,7)^«Cr 


-65 


3 


40Ca(a,7)44Ti 


-65 


27Al(a,n)30p 


-56 


4 


20Ne(a,7)24Mg 


-11 


30Si(a,n)33s 


-39 


5 


30Si(p,7)3ip 


-9.6 


i2C(a,7)^^0 


+19 


6 


36Ar(a,p)39K 


-7.5 


40Ca(a,7)44T- 


+15 


7 


27Al(a,p)30Si 


-4.0 


5«Ni(a,7)62Zn 


-8.7 


8 


33S(p,7)34ci 


+3.8 


27Al(p,7)28Si 


+6.0 


9 


i6O(a,7)20Ne 


-3.8 


24Mg(a,7)28Si 


+6.0 


10 


30Si(a,n)33s 


+3.5 


39K(a,p)42Ca 


+5.3 
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Table 8. Fits to thermally-averaged cross sections Na < av > (in cm^ s ^ mole ^) 

used in our calculations 



Reaction 


A 


B 


C 


D 


E 


F 


G 


*^Ti(a,p)''^V 


110.00000 


-3.84890 


-47.13600 


-91.31700 


7.65930 


-.54662 


30.27300 


Q(2a,7)^2Q 


CF88 














*5V(p,7)*«Cr 


90.47900 


.55899 


-37.67100 


-57.76700 


5.08410 


-.31734 


10.46700 


'*''Ca(Q,7)'*4Ti 


119.01000 


-4.29580 


105.93000 


-250.83000 


10.96000 


-.50058 


128.14999 


'5''Co(p,n)5^Ni 


13.38000 


-46.93500 


-3.84000 


9.20780 


-.27549 


.00430 


-3.93210 


^^Ar(a,p)3^K 


63.29300 


-14.78800 


-29.48900 


-42.76300 


5.14390 


-.42968 


11.34400 


''*Ti(a,7)4«Cr 


81.02500 


4.26030 


-245.58000 


137.92000 


-2.79730 


-.07990 


-98.51900 


i2C(a,7)i«0 


1.7x CF88 














"Ni(p,7)5»Cu 


109.94000 


-1.16540 


35.84500 


-155.78999 


9.22410 


-.50439 


64.77700 


5*Cu(p,7)5^Zn 


118.27000 


-1.34030 


45.34600 


-175.61000 


10.46300 


-.59085 


73.47000 


3^Ar(a,7)«Ca 


273.35001 


-5.81920 


235.08000 


-552.10999 


34.32800 


-2.01860 


245.53999 


^^Ti(p,7)^^V 


81.78900 


.07713 


-30.83200 


-55.66900 


5.66600 


-.43928 


10.56000 


5^Co(p,7)5**Ni 


70.29000 


-1.41470 


30.19200 


-105.66000 


5.11650 


-.27513 


50.64400 


5^Ni(n,7)5**Cu 


15.67600 


-.01193 


.42180 


-.30912 


.08421 


-.01823 


.01975 


^*Fe(a,n)5^Ni 


17.59100 


-67.55000 


1.68610 


-4.90450 


.88835 


-.05892 


1.50260 


''°Ca(Q,p)«Sc 


52.60400 


-40.91600 


-22.86600 


-28.98800 


1.59900 


-.06980 


10.16200 



Na < (Jv >= exp(A + B/Tg + C/tI''-" + DTg'^ + ETg FT^''' + GlnTg) 



1/3 



The rates for these reactions are from CF88 (CJaughlan & Fowler 198^) or arc from the formula 
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Fig. 1.— Final '^'^Ti mass fraction in a zone expanding with p oc as a function of initial 
density (g/cc) and temperature (K) for initial composition consisting of pure ^^Si. The contour 
levels show the logarithm of the final ^^Ti mass fractions relative to the final ^^Ti mass fraction 
in our reference calculation that began with initial temperature T9=5.5 and density p=10^ g/cc 
(+ sign). This figure shows that the initial temperature and density chosen for the survey in the 
rest of this paper (+ sign) produces a final ^^'pj j^^gg fraction that is roughly equal to the final 
^^Ti mass fraction over a quite large range of initial temperatures and densities. Most of the plane 
is within a factor 10. The dashed line shows the temporal evolution of temperature and density 
in the p oc expansion of our survey. The dash-dotted line shows the peak temperatures and 
densities of matter in the ejecta of the 20 Mq supernova model of Thielemann et al. (1990). The 
X sign locates the mass cut of the model at ~1.6 Mq radius in order to eject 0.07 Mq of ^^Ni. It 
is important to note that unlike the dashed curve, the dash-dotted curve is not a time history of a 
single zone. Rather, it shows the peak conditions of p and Tg reached by a continuum of mass zones 
in the explosion. The conditions chosen for our reference calculation are reasonably close to those 
achieved in supernova models. The ridge extending from the bottom left to the top middle of the 
figure comprises initial conditions for which the final ^^Ti mass fraction reaches its highest yield; 
this is the condition where the abundance of heavy nuclei within QSE is much the same as the 
abundance of heavy nuclei within NSE from the moment when quasi-equilibrium is first achieved 
until the freezeout. 
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Fig. 2. — The evolution of the mass fractions of some important nuclei in the nearly adiabatic 
expansion experiencing alpha-rich freezeout with standard nuclear reaction rates and initial 
temperature T9=5.5 and density p=10^ g/cc. 
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Fig. 3. — The ratio of the final Ti mass fraction to the "^^Ti mass fraction from the standard 
network if the reaction rates of ah charge-increasing and mass-increasing reactions of all isotopes 
of a single element are multiplied, element by element, by 0.01(a) or 100 (b). 
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Ar Isotopes 



Ti Isotopes 



31 32 33 34 35 36 37 38 39 40 41 42 43 44 41 42 43 44 45 46 47 48 49 50 51 

A A 



Ca Isotopes 




V Isotopes 



35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 

A 



42 43 44 45 46 47 48 49 50 51 52 53 54 
A 
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32 33 34 35 36 37 38 39 40 41 42 43 44 41 42 43 44 45 46 47 48 49 50 51 

A A 



Ca Isotopes 



V Isotopes 



35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 

A 



42 43 44 45 46 47 48 49 50 51 52 53 54 

A 



Fig. 4. — X{'^^Ti) /X{^'^Ti) standard ^^v isotopic reaction rates decreased by xO.Ol (a) and increased 
by xlOO (b). Reactions with ^Ti and ^^V produce the largest abundance sensitivity. 
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Fig. 5. — X{'^'^Ti)/X('^'^Ti)gia,ndard foi^ reaction rates decreased by xO.Ol (top) and increased by 
xlOO (bottom). Reactions with ^^Ti and ^^V produce the largest abundance sensitivity. 
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Multiplication Factor on *^V(p,7)*®Cr Reaction Rate 



Fig. 6. — Final X(^'^Ti) dependence on the reaction rate of ^^V(p,7)^^Cr. 
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Multiplication Factor on 3a Reaction Rate 
Fig. 7. — Final X("^^Ti) dependence on the reaction rate of a(2a). 
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Multiplication Factor on *'^i{a,y)*^Cr Reaction Rate 

Fig. 8. — Final X(^'^Ti) dependence on the reaction rate of ^^Ti(a,7)^^Cr. 




Multiplication Factor on **Ti(o(,p)*^V Reaction Rate 
Fig. 9. — Final X("^^Ti) dependence on the reaction rate of ^"^Ti(a,p)^^V. 
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Fig. 10. — Evolution of the Ti mass fraction as the temperature falls, with the standard triple- 
alpha rate (dashed curve) and also with the rate increased by xlOO (solid line). The dotted line 
shows the NSE mass fraction of ^^Ti during the expansion. 



-34- 




Fig. 11. — Evolution of the Ti mass fraction as the temperature falls, with the standard 
^^Ti(Q;,p)^^V rate (dashed curve) and also with the rate increased by xlOO (solid line). The 
dotted line shows the QSE mass fraction of ^^Ti during the expansion. 
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Fig. 12.— Evolution of the ''^Ti mass fraction as the temperature falls, with the standard 
^^V(p, 7)^*'Cr rate (dashed curve) and also with the rate increased by xlOO (solid line). The 
dotted line shows the mass fraction of "^"^Ti during the expansion if {p,-y) — {'y,p) equilibrium for 
the N = 22 isotones were able to persist. 
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Fig. 13. — Evolution of the (p, 7) — (7, J*) equilibrium for the N=22 isotones as the temperature 
falls, with the standard ^^V(p,7)^^Cr rate. The ^^Ti and ^^V pair (Z=22 and 23) increasingly 
depart from equilibrium with the ^^Cr and ^^Mn pair (Z=24 and 25) as the temperature declines. 
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Fig. 14. — The percentage changes of ^^Ti mass fraction as a function of rj when the standard 
^^V(p,7)^^Cr reaction rate is decreased by xO.Ol (fihed squares) and increased by xlOO (filled 
circles). The importance of '^^V(p,7)^^Cr reaction rate to the ^^Ti production in alpha-rich freezeout 
condition declines at r/ > 0.0004. The dotted line (open squares) shows the proton abundance at 
Tg=2 as a function of rj. The drop of the proton abundance at > 0.0004 causes the diminishing 
importance of ^^V(p,7)^^Cr reaction rate for increasing neutron richness. 



